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Abstract 

We perform a well defined derivative expansion to obtain the time dependent 
effective theory for a BCS superconductor at finite temperature, using an arbi- 
trary curve in the complex time plane. Our expansion is unique, being free of any 
problems related to the order in which analytical continuation and the derivative 
expansion are made, or equivalently to the order in which the zero momentum and 
zero energy limits are taken. In other words Landau damping terms are present but 
do not lead to the pathologies found in standard approaches. We discuss the reason 
for this, suggesting that these standard approaches may be internally inconsistent. 
The methods presented here are an attempt to overcome those inconsistencies, and 
finally, our results are compared to others in the literature. 



1 Introduction 

Effective actions have many applications in physics, since they simplify the description 
of highly complicated physical systems by integrating out all but the relevant degrees 
of freedom. As a result one is left with effective equations for the physically significant 
modes, with coefficients and couplings depending on the rest of the content of the theory. 
In the case of superconductors, static phenomena are well described by the effective 
Lagrangian of Ginzburg and Landau [1], and this can be derived from the microscopic 
BCS theory (see for example [2, 3] and references within). In relativistic particle physics, 
effective actions play an important role in inflationary theories of cosmology [4] , and can 
be used to describe low energy QCD and electro-weak phenomena [5] to name but a few. 
In the simplest case of fields which are constant in space and time, the effective action 
reduces to the effective potential multiplied by a space-time volume factor. 

However, to study the non-equilibrium dynamics say of a phase transition, be it in 
superconductors, inflation, or QCD, space and time dependent corrections to the effective 
potential are required. To obtain these, one often carries out a derivative expansion, the 
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lowest order of which gives the effective potential. Calculations of the derivative expansion 
have been successfully carried out at zero temperature T = where there is no Landau- 
damping, e.g. see [3, 6, 7] in the case of BCS superconductors, and references [5, 8, 9] 
for relativistic contexts. But at finite temperature, T > 0, the problems arising from 
Landau- damping have led to a failure of the derivative expansion both for BCS theory 
[10, 11, 12, 13] and in relativistic models [14, 15, 16, 17, 18, 19, 20, 21, 22]. 

So far, suggested solutions to this problem have been limited. One common approxi- 
mation which is often used to produce an effective theory is simply to replace the classical 
potential by the finite temperature effective potential, keeping the same classical kinetic 
terms. Indeed this approach is often used in simple models of inflation [4] and bubble 
nucleation in first order phase transitions (see for example [23]). It is very difficult to 
justify such methods, however, as the derivative expansions which should be used to ver- 
ify them fail at T > as commented above. In the case of BCS superconductors, Stoof 
[11] appreciated the problems but simply dropped the Landau damping terms. The same 
is true of [2]. Aitchison and Lee [13] try another approach. They include the effect of 
thermal damping, i.e. resum some effects from higher order terms into the propagators 
[13]. This solves the technical problem with the derivative expansion and is physically 
realistic, but experience suggests that it can be difficult to implement such resummations 
consistently [24, 25]. It is also technically very difficult to work with these methods. 
Many other authors have also tackled this question (see [3] and references within). 

The purpose of this paper is to study once more the question of derivative expansions 
at T > 0, and here it is done in the context of BCS theory. As mentioned above, this 
question has been the subject of numerous papers, not only in the case of superconduc- 
tors [3, 10, 11, 12, 13], but also in the apphcation of derivative expansions to particle 
physics and cosmology [14, 15, 16, 17, 18, 19, 20, 21, 22]. Despite this, we feel that this 
paper is justified as it raises some questions and potential problems with those previous 
calculations. Indeed, we will argue that we believe those calculations to be internally in- 
consistent and therefore incorrect. It is then argued that this inconsistency is responsible 
for the badly defined and non-unique T > derivative expansions which were previously 
obtained. 

The paper is set up in the following way. In section 2 the effective action for the pair 
field is derived from the BCS Lagrangian. This effective action can only be evaluated 
within an approximation scheme, whose first component is a weak field expansion (section 
3) . The second component is a derivative expansion (sections 4-6) . Section 4 discusses the 
reasons for which we believe the previous calculations of the derivative expansion to be 
internally inconsistent, whilst section 5 presents our own calculation scheme. Throughout 
we work at finite temperature and in arbitrary path ordered approaches to thermal field 
theory — that is, we do not work specifically in the imaginary- or real-time formalism, but 
with an arbitrary contour C in the complex time plane. As usual, the contour starts at 
ti and ends at ti — i(3 (where [3 = l/T), and our results are explicitly shown to be the 
same for any C. To carry out the calculation we work throughout in {t,p) space, and 
furthermore on its implementation in section 5 we see that Landau-damping no longer 
leads to problems with our derivative expansion which turns out to be well defined. The 
explicit results of our calculation are given in section 6. There we work for simplicity in 
the degenerate gas approximation [2] and present results for the effective potential as well 
as higher order corrections in time and space gradients. From those results the effective 
equations of motion for the Cooper pair field are seen to break time translation invariance, 
having an explicit dependence on the initial time t^ at which the system was set up. We 
believe this to be reasonable as by the very act of performing a derivative expansion one 



is dealing with time-dependent non- equilibrium systems. Our effective equations are also 
seen to contain both oscillatory and dissipative terms for all times t > ti. Finally, though 
we believe that the method presented here is an improvement on previous ones, it raises 
some confusing issues and may still not be entirely satisfactory. Such issues are discussed 
with the conclusions in section 7. 



2 The effective action 

In this section we derive the effective action for the auxiliary (pair) field A [2, 3]. We 
work in Minkowski space in (3+1) dimensions and use the notation x — {t,x) where, as 
usual in finite temperature calculations, the times can be complex. The starting point is 
the microscopic BCS Lagrangian 



(2.1) 



Here ips are anticommuting fermionic field operators, m is the electron mass and ^ is the 
chemical potential. The last term leads to the formation of bound states (Cooper pairs) 
for g > As usual, it is useful to work with the two component fields 



so that in equilibrium at finite temperature T — 1/ (5 the partition function is 
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Here d^x = dt (i^x where C is any contour in the imaginary time plane starting 
at ti and ending at ti — i[3 [26, 27]: here we do not restrict the calculation to any particular 
thermal field theory formahsm. In this path integral representation the *'s are now fields 
(and no longer operators) which are anti-periodic; 

*(t,x) = -*(t-i/3,x). (2.3) 

In the usual way, one can rewrite the quartic term in (2.1) through the identity 
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which on substitution into (2.2) and performing the Gaussian integral over the Grassman 
fields gives 



Here 



S,s[A*, A] = S,4A*, A]- - [ d^x A* A 

9 Jc 



where the effective action S'efr is given by 

5'eff[A*,A] = -iTrlnK-^ ^ -itr [ d'^x\nK-\x,x) 

Jc 



(2.5) 
(2.6) 

(2.7) 



and the matrix K ^ by 

Finally, in equation (2.7), 'tr' means the trace over matrix components. 

The purpose of the remainder of this paper is to calculate Ses- Such a calculation 
must clearly be done within an approximation scheme, and this consists first of a weak 
field expansion (section 3) and then a derivative expansion (sections 4 and 5). 



3 The Weak Field Expansion 

The first component of the approximation scheme is a weak field expansion. We work in 
the symmetry broken phase, so that < T < Tc where Tc is the critical temperature, 
and then expand A(i, x) in the effective action 5'eff[A*, A] by writing 



A(i,f) = (|Ao| + 0(t,f))e^^°. 
Here Aq = |Ao|e*^° is the x independent solution of the gap equation 
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One can absorb the phase by making a symmetry transformation so without loss of 
generality we will take Aq real, 6q = 0. 

In order to expand about Aq, first introduce the matrices $(a;) and G^^{x, y) defined 

by 



*(x) = 

G~\x,y) = 
and also the inverse matrix G(x, y) which satisfies 
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(3.3) 



j dS G-\x,y) G{y,z)^b\x-z)^. 
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Now, by definition of K ^ in (2.8), it follows that 

K-i(x,y) = G-\x,y)-r^{x)8\x-y) 

= y" d^z G-^(a;, z) {b^{x - z)n + G(z, y)^{x)\ . (3.5) 

Hence substitution of (3.5) into (2.7) followed by an expansion of the logarithm in powers 
of 0(x) gives 



5eff[A*, A] = S^^^W. 0] + S^^W. 0] + 



(3.6) 
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where 



Jc 

= -ill I d^xG{x,x)^{x) (3.7) 
Jc 

S^^WA] = '-It jj'x jj'yG{x,y)^{y)G{y,x)^{x). (3.8) 

Notice that whilst the first two terms in the expansion of S^q are local, the third one S),^ 
is non-local, as indeed are the higher order terms in 0. The purpose of the derivative 
expansion, which is the second ingredient in the approximation scheme used to calculate 
Ses, is to remove this non- locality. Here we focus on S^q , and the derivative expansion 
is done by Taylor expanding ^{y) in (3.8) as 

^{y) = ^{x) + {y- x)—^ + . . . . 

It is this derivative expansion which brings us into contact with the technical problem 
of ill-defined zero energy and momentum limits of thermal Green functions, and these 
in turn cause the inconsistencies usually found in such expansions at finite temperature 
[10, 11, 12, 13]. 



4 The derivative expansion 

Before carrying out the above mentioned derivative expansion on S^q, we try to sum- 
marise the different approaches to this expansion that have previously been used in the 
literature [2, 3, 10, 11, 13]. 

First consider the procedure carried out by authors who have worked in the imaginary 
time formalism of thermal field theory [2, 3, 10, 13]. Typically their calculations proceed 
through the following stages. Normally the first step is to Fourier transform equation 
(2.7) into momentum space to give 

S'S=-\^'Y. J (SsE J^,^(pMk-p,M-p)G(k). (4.1) 

ko Po 

A simple but important point to note here is that the energies, kQ,pQ, have been taken 
to be imaginary and discrete as is typical for systems in equilibrium. The second step 
is to expand expression (4.1) about p = in powers of po and p, and then evaluate the 
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sums and integrals over k. Finallv- the resulting expressions are converted back into real 
times (that is some analytic continuation is carried out) to give an expression for S'^g in 
terms of derivatives of ((){x). Notice that in the first step of these calculations — that is, in 
writing down equation (4.1) with the discrete momenta — both space and time translation 
invariance have been assumed, and furthermore the fields and propagators involved in 
the rest of the calculation are those of equilibrium systems. It follows that in the second 
step, discrete momenta (or equivalently equilibrium fields) have been used within the 
derivative expansion itself — below we will comment further on this point. 

The aim of this section is to discuss whether or not the operations outhned in the 
above paragraph are in fact consistent. 



To that effect consider now expression (3.8) in (t,x) space (rather than the Fourier 
transformed expression (4.1)). In particular, focus on the time component of (3.8) since 
it is this which changes between the T — and T > calculation, and also let x — (t, x) 
and y — (t',y). The derivative expansion of is obtained by writing 

^{t')^^{t) + {t'-t)^ + ... (4.2) 

Substitute (4.2) into (3.8) and consider the contribution to S'^g which arises say from 
the term containing the time derivative dt^{t). Clearly this term is not periodic in t' in 
the imaginary time direction since t' — t — i[3 ^ t' — t. Indeed this conclusion holds for 
every individual term in the derivative expansion of (2.7), apart from the ^'-independent 
contribution. 

Thus we arrive at the conclusion that each term in the derivative expansion of (3.8) is 
not periodic in imaginary time. (One should perhaps not be surprised at this — by the very 
essence of the derivative expansion one is dealing with slowly time varying and therefore 
non-equilibrium systems.) Hence it follows that at each order in the derivative expansion, 
taking a Fourier transform would not give rise to discrete energies, and hence it is not 
correct to use discrete energies (which correspond to periodic fields) if one were to work 
with Euclidean energies. However, we have stressed in the first paragraph of this section 
that authors who previously calculated derivative expansions in the ITF formalism did 
assume such discrete energies. We therefore beheve that these previous calculations are 
inconsistent. 

Before suggesting a method to overcome this inconsistency, note that the above dis- 
cussion has emphasised the imaginary time approach to thermal field theory whilst some 
authors have used a real time approach [11]. Recall though that these alternative real- 
time methods always give the same answer as the imaginary time calculations for any 
physical quantity [28, 29, 30]. In particular, this is true when one compares results for 
thermal effective actions, provided note is taken of various technical issues [19, 31, 32]. 
We now rephrase the above arguments in a way suitable for real-time methods. This is 
done by discussing the periodicity of the fields — a boundary condition which is common 
to all path-integral approaches to thermal field theory. The generahsation of expression 
(4.2) to infinite order is 

<^(t') = e^^"(*'-*)$(t") (4.3) 

t"=t 

where 

E" = -i A. 
dt" 

If the left hand side, ^*(t'), is periodic in imaginary time, then for equality the right hand 
side must also be periodic. This is assumed in previous calculations [11] and is clearly 
only true if e^^ = 1, a condition which was imposed in those calculations. However, as 
was noted above, each individual term in the derivative expansion (in other words the 
expansion of the exponential in (4.3)), is not periodic. Furthermore, since one always 
truncates the series and works to a finite order in derivatives, it follows that the approxi- 
mation to the field $(t) used in the expansion cannot be periodic. Thus, we believe that 
setting e^^ = 1 is not consistent with the derivative expansion. Again, we stress that 
since this was done in previous calculations, we believe these to be inconsistent. 

The approach to the derivative expansion presented in this paper therefore has one 
crucial difference with previous calculations. As shown above, each term in the derivative 
expansion is not periodic, and therefore our calculation does not assume periodicity. This 



is equivalent to saying that we allow e^^ 7^ 1, or equivalently that we do not work with 
discrete momenta for the field — as is typical of non- equilibrium systems. Given 
the lack of discrete momenta, it is difficult to work in (po,p) space and instead we work 
in {t^p) space. This is in fact a great simplification as it turns out that all the time 
integrations can be done exactly for any arbitrary thermal field theory formalism, and 
thus it avoids all the above mentioned issues. Furthermore, there are no problems with 
analytic continuation, which is another point we feel has been passed over rather rather 
hastily in the existing calculations. Finally, note that consistently with these discussions, 
we obtain the same results as in the existing hterature [10, 13] if we set e^^" = 1.^ 
Our convention for Fourier transforms is 



5 The calculation 

For the reasons discussed in the previous section, namely that each term in the derivative 
expansion leads to non-periodic functions, we work throughout in {t, k) space. The first 
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part of the calculation, section 5.1, consists of determining the propagator G(i, i', k) which 
appears in S).^ . The gap equation is then calculated in section 5.2 whilst the derivative 
expansion of S^^ is set up in section 5.3. There we will see that by working with the 
time variables, the integrations can indeed be done explicitly for an arbitrary thermal 
field theory formalism. At the same time this avoids the issues of analytic continuation 
which comphcate some approaches [28, 30, 33]. 



5.1 The fermion propagator 

As was discussed in section 4, the fermion propagator G(t, t', k) is in fact an equilibrium 
propagator since the Fermi fields are in equilibrium by equation (2.3). Its calculation 
therefore follows standard methods which we outline here; the details are given in ap- 
pendix A. 

First use (3.3) to define the matrix G~^{x) through 

G-\x)5\x-y) := G'^x^y) 

so that in (t, k) space 

G-\trk) = {'^'-^^^ \ = {'^'z'' -J^i ) (5.1) 

where eg is the energy measured from the Fermi surface; 



^Regarding the discussions of (anti-)periodicity made above, it should be noted that G{t, k) is the 
propagator for the Fermi fields, which by definition in equation (2.3) are anti-periodic in imaginary time. 

Hence G(t, k) is anti-periodic; G{t, k) = —G{t—i(3, k). Regarding fc), or equivalently the field A(t, fc), 
we have shown that for the derivative expansion to be consistent at any finite order, A{t, k) cannot be 
periodic. That this is the case could be achieved by adding a small time dependent perturbation to A, 
and wc arc then watching the system respond to this. In addition, it should be noted that in equation 
(2.4) which introduces the A{t, k) field, there is nothing to constrain A to be periodic in the first place. 
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As usual in the BCS limit in which we are working here, the Fermi wavenumber is defined 
through the chemical potential fihj kp = 2m/i. 
Now, from equation (3.4), G{t,t',k) satisfies 

G-\t,k)G{t,t' ,k) ^ 5{t - t')1L. 

In the usual way one can solve this equation for G{t, t', k) by diagonalising the operator 
G~^{t,k) which given in (5.1) (see appendix A). The result is 

Qu t' k)=( <F_{t,t',k)+vlF+{t,t\k) u+u.{F+{t,t',k)-F.{t,t',k)) 
^' ' ^ yu+u^{F+{t,t',k)-F_{t,t',k)) ulF+{t,t',k)+u'iF_{t,t',k) 

' (5.2) 

Here u± are functions of k (this dependence has not been written in (5.2) for clarity), 
and are given by 









i±ii 







u±{k) - + 
where is the dispersion relation 

It is important to recognise that u+{k) and U-{k) arc just the usual functions often 
called u{k) and v{k) respectively which, in the operator approach, are the coefficients 
of the Bogohubov transformation which diagonahses the Lagrangian. We have chosen 
this different notation as later (in section 5) it will enable us to write our results in a 
significantly more compact form. 

As a result of the diagonalisation procedure (appendix A), the functions F±{t,t',k) 
which appear in (5.2) are seen to satisfy the Feynman equation 

{i^t±E^)F±{t,t',k)^s{t-t')^L. 

Hence their solutions are known, and for our purposes they are best given in a Mills 
representation [34, 32]; 

F±{t, t\ k) ^ -ij ^e-*(*-*') {ec{t, t') - N{ko)) p±{ko, k). (5.3) 

Here N{ko) is the Fermi-Dirac distribution: 

the 9c is the usual 9 function for the contour C [26, 27], and the spectral functions p are 
given by 

p±{ko,k)^2nS{ko±E^). (5.4) 
Finally we note that F+{t, t',k)^ F_{t\ t,k). 
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5.2 The Linear terms and the Gap Equation 

Before turning to the calculation of in (3.8), we consider briefly the gap equation (3.2) 
or in other words the terms linear in in S^s- As usual the gap equation is a necessary 
tool in the evaluation of S^^ as it is used to remove UV divergences (see below). 
Equation (3.7) gives S^^ to be 

SiS[A*,A]^-i f d''x{Gum*{x)+G2im{x)) 
Jc 

where from (5.2) 



6-12(0) :=Gi2(a:,x) 



(27r)3 2E^, 
d^k Ao 



{F+{t,t,k)-F_{t,t,k)) 
l-2N{E^)). 



Thus on substitution into (3.2) one obtains 

d^k 1 



(5.5) 



which agrees, for example, with Stoof [11]. Notice that the first term on the left hand 
side of (5.5) is linearly divergent since E^ oc k'^ for large k. 



5.3 Quadratic Terms 

Now consider the important quadratic term of the effective action (2.7). Substitute 
the propagator (5.2) into (3.8) Fourier transformed into (t, k) space. On taking the trace 
in (3.8) one obtains four terms. Two of these are proportional to (f)*(t,p)(f)(t' , —p) and 
(}){t,p)(})*{t' , —p) and, as is common in the literature, we call these the diagonal terms. 
The non-diagonal terms are proportional to (f)*{t,'p)(j)*{t' , —p) and (l){t,p)(j){t' , —p) and in 
fact the coefficients of these two non-diagonal terms will be seen to be identical. Thus 
we write 



S. 



. (J = Si^ [Da] + [Db] + Si^ [ND] (5.6) 

where the D and ND stand for diagonal and non-diagonal as explained above. The 
diagonal terms are given by 



dt / dt' 
2 Jc Jc 



d^p 



d^k 



0(t,p)0*(t',-p)x 



X 



X 



X 



and 
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(27r)3 J (27r)3 

ul{k - p/2)F+{t, t', k - p/2) +u^_{k- p/2)F_(i, t', k - p/2) 
uiik + p/2)F_ {t', t,k + p/2) +u^_{k+ p/2)F+{t', t,k+ p/2) (5.7) 



ul{k + p/2)F+{t', t,k+ p/2) +u'i{k+ p/2)F_ {t', t,k + p/2) 
ul{k - p/2)F_{t, t', k - p/2) +ul{k- p/2)F+{t, t', k - p/2) .(5. 



Observe that interchange of t and t' in (5.7) followed by the relabelling p —p yields 



S'gg [Db] = Sgg[Da]. However, consider now any given term in the derivative expansion 
of ^£^[Da], which is obtained by writing -p) = (j)*{t, -p) + {f - t)dt(p*{t, -p) + . . . 
in (5.7). It is straightforward to see that the changes t' t and p — > — p for that term 
do not yield the corresponding term of [Db] (the only exception is the t' independent 
contribution). That is, whilst this change of variables will yield yS^g [Db] = Sgg[Da] to 
infinite order in the derivative expansion, this is not true to any finite order. Since the 
derivative expansion is always truncated, it is therefore necessary to work with S'^g [Db] 

('2) 

and S),^ [Da] separately. Below we will indeed verify that the results obtained in this way 
differ from those of a derivative expansion on 25*^0 [Db] (see equation (6.4)). 
The non-diagonal term in (5.6) is given by 



^(2) 



I dt I dt' I 

Jc Jc J 
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d^k 



{<f)*{t,p)<p*{t', -p) + <p{t,mt', -P)} 
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X u+{k- p/2)u- {k - p/2)u+{k + p/2)u- {k + p/2) 



X 



X 



F+{t',t,k+p/2) - F_{t',t,k+p/2) 
F+ {t, t',k- p/2) - F_ {t, t',k- p/2) 



(5.9) 



5.3.1 Derivative expansion of 5'eff[Da] 



As a specific example we now carry out a derivative expansion on 5*65 [Da] given in (5.7). 
The computation of 5'eff[Db] and S'eflf[ND] will follow very similar steps. 
First Taylor expand 0(i', —p) about t' — t; 



ct>{t\ -p) 



<f>{t, -p) + {f - t) 



d(l)it, -p) 



dt 



2! 



it' - tf 



d\t,-p)ct> 



dt'' 



t"=t 



+ ... (5.10) 
(5.11) 



where E" — —idf . It is now possible to proceed in two different but equivalent ways. The 

first is to consider a given term in (5.10) (for example one proportional to 9^0), substitute 
it into (5.7) and then calculate its contribution to the effective action. Equivalently one 
can substitute (5.11) directly into (5.7), and then obtain the contribution from by 
simply calculating the term proportional to in the expansion of [Da] in powers 
of E". We follow this second more powerful method here. 
Therefore, substitute (5.11) into (5.7) to obtain 



dt 



c 



d^p 
(2^ 



(l>*{t,p)BDa{t,ti,t",p)(l>{t",-p} 



where Boa is given by 

i f d^k f 
S..(M.t",p-):=-y^ J^dt' 

X \ul{k - p/2)ul{k + p/2)F+(i, t\ k - p/2)F+{t', t,k+ p/2) 

+ ul{k- p/2)ul{k + p/2)F4t, t', k - p/2)F4t', t,k + p/2) 
+ ul{k - p/2)ul{k + p/2)F+{t, t', k - p/2)F_{t\ t,k+ p/2) 



+ u^_{k- p/2)ut{k + p/2)F4t, t', k - p/2)F+{t', t,k + p/2) 
- U-^^ H u'.Jk-mu-Jk + mZsosAt,U,E'',k,p) (5.12) 



and 



Z,,s^it,ti,E",k,p):= I dt' F_so{t,t',k-p/2)Fs,{t',t, 

Jc 



k+p/2)e 



iE"{t'-t) 



(5.13) 



In traditional work [10, 13], one finds that the retarded thermal Green function [28, 29] 
appears instead of Boa- These are the same only in the limit e^^ = 1 in accordance 
to the discussion in section 4. Below we will see that the so called Landau damping 
contributions come from the terms in (5.12) with Sq = — Si, and furthermore that they 
do not give rise to any problems with the derivative expansion. 

The only dependence on the variable t' in (5.7) is now in Z, and we show that this 
may be evaluated for any contour C. Thus our method includes both imaginary- and 
real-time approaches to thermal field theory. Substitute expressions (5.3) for F± into 
(5.13) to obtain 



Z.. 



)S1 



dt' I ^ [ ^ p^so{ko,k-p/2)PsAki,k + p/2) x 



dko 
2^ 



27T J 2t\ 

N{ 
dki 

2^ 



c 



p-soiko, k - p/2)ps,{ki, k + p/2)Jc{ko, ki, U - t, E")(5.U) 



Observe that all the time dependence is in the function Jc{ko,ki,ti — t,E") which, on 
letting A — E" + ko — ki, is given by 

Jc{ko,k^,ti-t,E") = [ dt' [ec{t',t)-N{ki)][ec{t,t')-N{ko)]e'^^''-''^ 

Jc 

= / dt'N{ko)N{ki)e'^^''-'^ - N{ki) f dt'ec{t,ty^^''-'^ 
Jc Jc 

dt'9c{t',t)9c{t,t')e'^^''-'\ 



N(ko) f dt'ec(t',t) 
Jc 



C 



Consider now an arbitrary contour C starting at ti and ending at ti — ij3. The integrals 
in (5.15) may be calculated explicitly for that contour C: using the definition of 9c{t',t), 
the final term contributing to Jc{ko, ki, ti — t, E") vanishes; the only contribution to the 
second term is for ti < t' < t; and the contribution to the third term is for i < f < ti — i(5. 
Thus the second term, for example, is given by 



N{ki) I dt'ec{t,t')e 
Ic 



'\JA{t'-t) 



N{ki)e-'^' I dt' e 
It 



I JAt' 



After similar manipulations for the other terms, we obtain 



Jc(ko,ki,5,E")^ 



N{ko)N{ki) 

iA 



^Pk, ^^PA g/Jfeo (^l _ ^iA8\^ _ ^PE"y^ (5.15) 



where 6 = ti — t. This step is perhaps the most important one of our method — by working 
with the time variable t itself rather than Fourier transforming, we have been able to do 
the integrals explicitly and without using the explicit form for the spectral functions p. 
Substituting (5.15) back into (5.14) and using the definitions of p in (5.4) gives 

Zsosi {t, ti, E", k,p) — 

1 1 1 

(exp{/3so^fc_p72} + 1) (exp{/9siE,,^./2} + 1)^ 

X { (1 - e^^) + eMPisoE^,-,/, + s,E^,^.^,} (e^^^ - l) (l - e^^") } (5.16) 

where we have used the shorthand notation 

A = E" + soE^,_.^, + s^EJ;^./,. (5.17) 

Before going any further we make the following comments. First, observe that if we 
set e^^ = 1 then our expression (5.16) reduces to the corresponding equation (number 
(10)) of [13]. This is consistent with the discussion of section 4. 

Secondly, note that the time independent contributions to the effective action have 
E" — and hence e^^" — 1. Thus all our time independent results will be identical to 
those of [11, 13]. 

Thirdly, equation (5.16) has an explicit dependence on the initial time ti. Time 
translation invariance has therefore been lost as expected for such a time dependent non- 
equilibrium system. Furthermore, an important point to observe is that even a.t t = ti 
{S — 0), our expression (5.16) differs similar terms in [10, 13] by a factor of e^^" . This 
point will be discussed further below. 

Finally, consider the terms with sq = —si in Boa- From equations (5.12) and (5.13) 
and the expression for Z in (5.16), these terms correspond to Landau damping terms in 
the original formalism (when e^^ = 1). The reason is that in this case, the denominators 
contain E" ± {E^_^^2 ~ which vanishes when p — and E" — thus leading to 

a potential divergence. It is this divergence in the bubble diagram which is sometimes 
ignored [11], and was the cause of the badly defined derivative expansion in [13]. However, 
observe that in our case we have eJ-^^ 7^ 1, so that thep = and E" = limits are actually 
well defined: the zero in the denominator B^a when sq = —Si is exactly cancelled by a 
zero in the numerator of this expression. Thus here we are dealing with Landau damping 
terms that have a well defined derivative expansion. Similar comments will apply for the 
other diagonal term as well as the non-diagonal terms which we now calculate.^ 

5.3.2 Second diagonal term 'S'g^^[Db] 

Calculation of S^q [Db] follows a very similar route to that above. Indeed, the function 

Z introduced a])o\'e also a.])])eais in [Db] as can be seen by comparing the temporal 

^Interestingly, the cuts at E" = ±2Ao familiar from Green functions commonly encountered in field 
theory are also not present in the S's when sq = Si. However, one must be very precise when making 
the link between physics and the analytic properties of particular Green functions. For instance it is 
the theory of linear response which links cuts at E" = ±2Ao in retarded Green functions to physical 
decay processes at those energy scales [27, 29] . Here we are not using linear response theory but rather 
looking at systems which are slowly varying in time. It is therefore not clear that B should have the 
same analytic structure as the retarded Green functions of linear response [27, 29] as they found to be 
relevant to different physical problems. 

1 o 



dependence of 6*^5 [Da] in (5.7) with that of S'^g [Db] in (5.8). This leads to a switch from 
U- to u+ functions and one obtains 

si^[Dh] = r{t,p)Bnb{t,U,t",mt",-p} 

where the bubble diagram BDb{t,ti,t" ,p) contains contributions from all the four terms 
in (5.8): 

^^^=2/(2^ J] <(fc-p72)<(A^+p72)^3o... (5.18) 
5.3.3 Non-diagonal term -S^ff^[ND] 

The non-diagonal terms (5.9) again follow a similar pattern. Here, since 



u+{k-p/2)u+{k-p/2)u^{k-rp/2)u^{k-rm) = ' °' 

it follows that 

Jc J l^^j t"=t 
where the bubble diagram B]\fj:){t,ti,t" ,p) is 



B 



'ND 



\ I AT? ^^"l XI (-soSi)^.o.i- (5.19) 

2 J (27r)^ ^E^^p/2Ek-p/2 ,„^± 



5.3.4 Total quadratic term and derivative expansion 

To conclude, we have 

where 



t"=ti 



with the S's given in (5.12), (5.18) and (5.19). 

It is now straightforward to expand this effective action in powers of d'^cj) and 
where n > 0; it is done by expanding the B's above about E" — and p = 0. We now 
carry out this derivative expansion. 

6 The derivative expansion in the degenerate gas ap- 
proximation 

Though it is not difficult to expand the expressions given above in E" and p to obtain the 
derivative expansion and hence effective equations of motion, the resulting expressions 
are however very complicated. To simplify them we therefore work in the degenerate gas 



approximation [2] for which T < Tc, and furthermore the gap Aq is much bigger than 
the temperature. 

Before clarifying further this approximation, note that our expansion also imposes 
P\dtln{(f))\ <^ 1 since the factors of exp{PE"} = ex.p{—ipdt"} must be expanded. It is 
fundamental to remember that by their very definition, derivative expansions are valid 
only for fields varying slowly in space and time, and we should not be surprised that 
our expressions differ from those used in linear response theory. The latter uses retarded 
Green functions which describe the response to a delta-function impulse, a very different 
physical problem. Also with <S T for the expansion to be vahd, clearly taking 

the zero temperature T — > limit makes no sense in our expressions. 

To summarise, the results presented below are valid in the regime 

and the gap Aq will be used as the primary mass scale for the problem. Thus we introduce 
the following dimensionless parameters 

y = ±, .^ = (1 + ,^) = |!. 

To exploit the degenerate gas approximation, observe that the integrals are sharply 
peaked around the Fermi surface k = kj. This is not true if there are UV divergent 
terms but we shall see that in the effective potential, these term cancel through the gap 
equation. The Fermi-Dirac function Nf then appears as 

7V/ = [exp{/3Ao4 + l]-\ 

and is even more sharply peaked than other parts of the expression. Thus we perform 
the following simplifications :- 

• Isolated factors of k not inside an e^^ or Ej^ factor can be replaced by a constant kf. 
Valid since k'^ = kj{l + y'^ /Y^) where Y := h/Aq > 1. 

• In terms with at least one factor of Nf, ~ 0, E'g ~ Aq, i.e. y ~ 0, ^ ~ 1. Valid 

since Aq 3> T. 

The results then have a common factor of p 

mkf 



For example, the gap equation in this notation is 

The expressions will also depend on the Minkowskii time t elapsed since the initial con- 
ditions were set at time tj. 
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6.1 Results 



The notation in this section is as follows. Prom equations (2.6) and (3.6), the effective 
action is given by 



S. 



eflf 



+ 



-itr J d^x lnG-^(x,x) -^Jd^x 



(6.1) 



where (f) — (f){t, x), (f)* — (t){t, x) and 7's are of the form 



a,b ^ ^ 



(6.2) 



The zeroth order effective potential terms, 7^0 easily obtained by first setting 
E" = = p in (5.13) and then using the gap equation to remove the UV divergences. In 
the degenerate gas approximation we have 



7o,o 



Db _ ND 

/o,o ~ /o,o 



dy 



2Nf - 2{Aof3)Nf{l - Nf) 



(6.3) 



The first order terms are slightly more complicated in part and are given by 

■ 0, 



7o,i — 


^Db _ 

7o,i — 


^Da _ 


P/3 f 


7i,o — 






-N/ (< 


^Db _ 




7i,o — 


16 J 




-N/ (, 


^ND _ 


P(3 [ 


7i,o — 


16 J 



^ND 

7o,i 



dy 



{z-yf 



-2Nf e 



,-2i{t-ti)Aoz 



^2 ^^2i{t-U)A^z _ g-2i(t-tOAoz^ ^ 2(Ao/3)A^/(l -Nf) [1-2 



i{t -ti) 



[z + yf 



- 2Nf e 



-2i(t-ti)Aoz 



-2i{t-ti)Aoz^ 



+2(Ao/3)iVj(l - TV^) 1 - 2 



(6.4) 



Notice that the only difference between these two diagonal contributions 7^^q and 7]^q are 
the terms odd in y (i.e. e^). One should also note that for t > ti, the 71^0 terms contain 
both real and imaginary components. In the equations of motion, which follow directly 
from (6.1), these correspond to oscillatory and dissipative terms respectively. The second 
order terms may be found in appendix B. 

The above results may be rewritten in terms of the real Higgs 4>h and Goldstone fields 
(f)g where 

0= -^{<Ph + i<Pg)- 



1 n 



One obtains 



'eff 



(x, x) — - I d^x An 



+ 



-itr j d^x lnG-^( 



The coefficients 7^* and 7^ are defined in the derivative expansion as in 6.2 and their 
exphcit forms are 



= / dy [^e-2^(*-**)^°^ + 2{Aof3)Nf{l - Nf) 

8 J z 





= 270^0" 




= 


Tfo 


8 . 


lo,i 


= 




= ^0,1 = 


11,0 


8 . 


7o,o 


- _ 

— 10,1 - 


7i,o 


-'41 



dy 



- - 2iV; - 2(Ao/5)iV;(l - Nf) 



iit - tj 



^ j dy Q - 2NA e-2^(*-*^)^°^ - 7V/ (e 



^2i(t-ti)Aoz _ g-2i(t-ti)Aoz> 



7 Discussion and Conclusions 



(6.5) 



Several comments can be made about these resuhs. First, for any static field configura- 
tions E" = 0, and hence exp{/3£^"} = 1 as assumed in all the previous calculations (see 
the discussion in section 4). Thus for static configurations we recover the static limit of 
retarded equilibrium Green function results used by others, so that our effective action is 
identical to standard calculations in the static limit. In particular, our lowest order term 
in the derivative expansion, the effective potential, is identical to the one which would be 
calculated using pure equilibrium methods. The differences are in the non-static terms: 
we have a unique derivative expansion which does not depend on the order in which the 
spatial and temporal derivative expansions are made. Thus only by using our expansion 
can one justify the use of the traditional effective potential as the potential for some 
effective theory, since here we have a fully defined and controlled expansion. Following 
on from this, since the physical masses for the bosonic field A are encoded by the second 
derivative of the effective potential, we expect to reproduce Goldstone's theorem. We 
indeed find that the mass term for the Goldstone mode, 7^q, is zero. 

As discussed in section 4, our calculation takes into account the fact that time trans- 
lation invariance is broken by the non-equilibrium derivative expansion. This is refiected 
in the final effective action obtained in (6.1) where the coefficients depend explicitly on 
the initial time ti at which the system was perturbed. 

Though we believe that the methods presented here are more consistent than previous 
methods used to calculate the derivative expansion, there are still a number of confusing 
aspects in this work. One of these (which we believe also exists in the previous calcula- 
tions, though it is not explicitly commented upon there) is the question of the difference 



between [Da] and [Db] of (5.7) and (5.8). We believe to have argued correctly in 



:.(2) 



section 5 that the change of variables t ^ t' and p —p should not be carried out before 
the derivative expansion. However, this in turn is responsible for the coupling between 
the Higgs and Goldstone modes in (6.5). Another interesting aspect of our results is 
that apart from breaking time-translation invariance (which we do believe to be correct), 
observe that another symmetry has been broken:- since (71^0) 7^ 7i|o then to this order 
in the derivative expansion S^cfr is not invariant under 0^0*. However, one fully expects 
this to be broken if there is a non-zero charge in the system. This corresponds to having 
an imaginary initial time derivative for the scalar field (j). Finally, note that the first 
terms of and 7^5 are actually divergent at t — ti. This infinity is exactly the one 
which arises in (5.13) if one were to set t = ti and p = 0. It should be noted that this 
divergence is not present in the usual calculations. As commented on just after equation 
(5.17) it stems from the fact that aX t = ti our expressions differ from similar terms in 
[11, 13] by the factor of e^^" ^ 1. In this non-relativistic case such infinities are not too 
worrying one always has a UV cutoff in mind, but the same is not true in relativistic 
cases where they would have to be renormalised in some way [35]. 

To conclude, we have argued that the traditional approach to producing derivative 
expansions of thermal effective actions are inconsistent since they try to use equilibrium 
fields in constructing the derivative expansion. We have argued that this is technically not 
correct, and that it is incompatible with the physical problem being described — that of 
slowly varying fields. We have shown that letting the fields be out of equilibrium leads to 
a well behaved derivative expansion. It means that the Green functions at the heart of our 
calculations, our i3's, are not the retarded thermal Green functions normally encountered 
in such work. Indeed, as linear response analysis shows [27], those are relevant to a very 
different type of physical problem where the system is hit with a sudden change rather 
than responding to a slow variation. Finally, we have used these ideas to calculate the 
effective action for a BCS superconductor. 
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A Calculation of the fermionic propagator 

Calculation of the fermion propagator in (t, k) variables proceeds as follows. Throughout 
it is vital to keep track of the delta functions and boundary conditions so that the solution 
is valid for any path-ordered formalism of thermal field theory. 
As in section 5, define the matrix G~^(a;) through 

G-\x)5\x-y):^G-\x,y) (A.l) 

— * 

so that in (i, k) space 
G-^{tS) = 




Ao \ _ ( idt- Ao 

2m 



(A.2) 



Here e^^^-^i, 



( -^k ^0 \ _ Y - cos(2^,^) sin(2%) \ 



and Ej^ is the dispersion relationship = \Aq\ + ei. The angle 9^^ is defined through 
cos(2^,^) = -|, sin(2^g) = ^ = -2K(^)|K(^)| (A.4) 



where, as in section 5, 

u±{k) 



so that cos(^g) = u+{k) and — sin(6'^) = u_{k). 
Now, from equation (3.4), 











i±ii 











(A.5) 



{idt-S. + A^)G{t,t',k) = (A.6) 
The diagonahsed form of Aj^ is 

A,— B,^ . A,^ . ^ (A.7) 
where Aj^ contains the eigen-values of Aj^; 

= ( t 4 ^ • '^-'^ 

and Br its eigen- vectors; 



_, -sin(%) cos(%) \_( u_{k) u+{k) , , . 



cos(%) -sin(%) ; \^ u+{k) u_{k) 
Now equation (A.6) is equivalent to 

{zdtBz' + BZ'A^) G{t, t', k)B^, = 6{t - t')1L (A.IO) 
so that on defining the matrix 

F(t, t', k) = B-iG(t, t\ k)Bf: (A.U) 

one obtains 

{idt-K + Ag) F(t, t\ k) = 5(t - t')-S.. (A.12) 
Thus satisfies the usual Feynman equation and is given by 



where in the Mills representation [34] 

F±{t,t',k) = -iJ^e-''°^'-'">{9cit,t')-N{ko))p±{ko,k) 

so that F_{t,t',k) = F+{t',t,k). Here N{ko) is the Fermi-Dirac distribution: 

^(^o) = (A.14) 

the 9c is the usual 9 function for the contour C [26, 27], and the spectral functions p± 
are given by 

p±{ko,k) = 27r5{ko±Ej:). (A.15) 

Finally, these results may be used to compute the explicit form of G{t,t',k): inverting 
equation (A. 11), using (A. 9) and (A. 13) gives, after some manipulation, 

C(ff'k)=( ulF_{t,t',k)+ulF+{t,t',k) u+u4F+{t,t',k)-F_{t,t',k)) 
^ ' ' ^ \^ u+u4F+(t, t', k) - F4t, t', k)) ulF+(t, t', k) + ulF_(t, t', k) 

— * — # 

where the k dependence of u±{k) has been omitted but is understood. 

B Second Order Terms 

Here we write out the results for the second order terms. 

7?o^ = / dy + i^(/3Ao)e-^(-*')^-[(/3Ao - 2Aoi(t - U))z - 1] 

+Nf[-2 + {A/3){AoPf] - Ar/e^^(*-*^)^°^(/5Ao)[l + /?Ao - 2Ao^(t - U)] 

-4Nf{l - Nf){f3Ao){Aot{t - U))[(3Ao - Aot{t - U)] 

+ (2 - Ar^)Ar^e-2^(*-*^)^°^(/?Ao)[l - {/3Ao - 2Aot{t - U))]} 



(A.16) 



Db _ Dal 
72,0 — 72,0 1 



IQAl 



I dy + ^^±^(/3Ao)e-^^(*-*^)^°1(/3Ao - 2Aoi{t - U))z - 1] 

+Nf[-2 + (4/3)(Ao/3)^] - Ar/e2^(*-*^)^°^(/5Ao)[l + /?Ao - 2Aoz(t - U)] 
-4iV^(l - Nf){pA^){A^i{t - t,))[/5Ao - A^i{t - U)] 
+Nf{2 - iVy)e-2*(*-*')^°^(/3Ao)[l - (/?Ao - 2Ao^(t - U))]] 

+Nf[2 + (4/3)(Ao/3)='] + iV/e2^(*-**)^°^(/?Ao)[l + /3Ao - 2Aoi{t - U)] 
-ANf{l - Nf){l3Ao){Aoi{t - UWAo - Aoi{t - t,)] 
-Nf{2 - Ar^)e-2^(*-*^)^°^(/?Ao)[l - (/3Ao - 2Aoi(t - U))]} 



Da _ Db 
70,2 — 70,2 



48mAg 



[-IQz^ + 10) + — y(-6z^ - 22>z^ + 20) 
A* 



+ANf [3 + 3(/?Ao) - (/3Ao)'] + 127V/(/3Ao) (-1 + (/3Ao))} 



-4z' + 10) + + 20) + 5 f^)' + - 2) 

A* V A* 



-47V^ [3 + 3(/3Ao) + (/3Ao)^] + 127V/(/3Ao) (1 + (/3Ao))} . 
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